Formation of long-lived, scarlike modes near avoided resonance crossings in optical microcavities 
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We study the formation of long-lived states near avoided resonance crossings in open systems. For three 
different optical microcavities (rectangle, ellipse, and semi-stadium) we provide numerical evidence that these 
states are localized along periodic rays, resembling scarred states in closed systems. Our results shed light on 
the morphology of long-lived states in open mesoscopic systems. 
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Understanding the properties of long-lived, quasi-bound 
states in open mesoscopic systems is of central importance for 
many research subjects, e.g. semiconductor ballistic quantum 
dots m 0. HI 01 J photoionization of Rydberg atoms |5], mi- 
crowave sy stents I aBl ; quantum chaos t8J, and optical micro- 
cavities lliHIIMl lilHQ. In several of these studies 
the long-lived states are scarred. The original scarring phe- 
nomenon has been discovered for closed systems in the field 
of quantum chaos lfT6ll . It refers to the existence of a small 
fraction of eigenstates with strong localization along unsta- 
ble periodic orbits of the underlying classical system. In open 
systems, however, scarred states seem to be the rule rather 
than the exception. The nature of the mechanism behind this 
scarlike phenomenon is not yet understood. 

Avoided level crossings in closed or conservative systems 
are discussed in textbooks on quantum mechanics. They occur 
when the curves of two energy eigenvalues, as function of a 
real parameter A, come near to crossing but then repel each 
other [17J. This behaviour can be understood in terms of a 
2x2 Hamiltonian matrix 

For a closed system this matrix is Hermitian, thus the energies 
Ej are real and the complex off-diagonal elements are related 
hy W ^V* . The eigenvalues of the coupled system, 

E1 + E2 I (El - E2y 

E±{A)^ ' ^ ' ±^ '^ ' ^ " +VW , (2) 

differ from the energies of the uncoupled system Ej only in a 
narrow parameter region where the detuning from resonance, 
\Ei{A) — i?2(A)|, is smaller or of the size of the coupling 
strength y/VW. The parameter dependence of V and W can 
often be safely ignored. 

The matrix ([T]l also captures features of avoided resonance 
crossings (ARCs) in open or dissipative systems if one allows 
for complex-valued energies Ej. The imaginary part deter- 
mines the lifetime Tj cx 1 /lm{Ej ) of the quasi-bound state far 
away from the avoided crossing i?2p S> l^VF, where the 
off-diagonal coupling can be neglected. Keeping the restric- 
tion W = V* allows for two different kinds of ARCs l isll . 
For 2\V\ > |Im(£^i) — Im(£'2)|, there is an avoided cross- 
ing in the real part of the energy and a crossing in the imag- 
inary part. At resonance Re(i?i) = Re (£'2) the eigenvec- 
tors of the matrix ([T]i are symmetric and antisymmetric su- 
perpositions of the eigenvectors of the uncoupled system. If 



one of the latter corresponds to a localized state then such an 
ARC leads to delocalization and lifetime shortening 1 19]. For 
2\V\ < \lm{Ei) — Im(i?2)|, there is a crossing in the real part 
and an avoided crossing in the imaginary part. This kind of 
ARC has been exploited to design optical microcavities with 
unidirectional light emission from long-lived states ll20ll . 

The case W — V* is called internal coupling since the only 
difference to the Hermitian coupling of two states in a closed 
system is that each state is individually coupled to the con- 
tinuum. The latter is described by the imaginary part of the 
diagonal elements of matrix The fully nonhermitian case 
W ^ V* is more general; it permits an external coupling of 
the states via the continuum. Figure [T]illustrates that the real 
part undergoes an avoided level crossing as in the case of a 
closed system. The important feature is that one of the states 
has a considerably increased lifetime. The constraint of the 
conservation of the trace of the matrix in Eq. ([T]i simultane- 
ously generates a state with short lifetime. The formation of 
fast and slowly deca ying st ates is known as resonance trap- 
ping, see, e.g., Refs. [|2ll 12211 . 

The aim of this letter is to show that ARCs due to external 
coupling can have a strong impact on the localization proper- 
ties of long-lived states in open systems. The symmetric or an- 
tisymmetric superpositions are more localized in real or phase 
space than the original states, so that important decay chan- 
nels are blocked. A surprising finding is that these states can 
resemble scarred states which helps to explain the frequently 
observed scarring in open mesoscopic systems. 

We examine optical microcavities, where the optical modes 
and their frequencies play the role of states and their energies. 
Light confinement in microcavities has attracted considerable 
interest in recent years due to the huge potential for various re- 
search fields and applications, for a review see Ref. ]23]. For 
most applications, like low-threshold lasing, long-lived modes 
are required. We consider quasi-two-dimensional dielectric 
cavities with rectangular, elliptical and stadium-shaped cross 
section. We first focus on rectangles because of the following 
convenient properties: (i) the modes not close to an ARC can 
be computed analytically to a good approximation, including 
mode pattern and complex frequency; (ii) the internal ray dy- 
namics is trivial, so localization effects related to chaotic ray 
dynamics ll24ll can be ruled out. Rectangular and square mi- 
crocavities have already been studied both experimentally and 
theoretically ]25, 26, 27]. However, ARCs in these systems 
have not been addressed so far 

We fix one side length to R = 2/im and vary the aspect 
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FIG. 1: Avoided resonance crossing in the case of the matrix (TJ with 
VW = 10.000225, El = -iO. 015, and £'2 = A-iO.Ol. Real (top) 
and imaginary (bottom) part of the energy vs. parameter A. 



ratio < £ < 1. We choose the effective index of refrac- 
tion to be n = 3.3 inside and n = 1 outside the dielectric for 
the transverse electric (TE) polarization with magnetic field 
H perpendicular to the cavity plane. Maxwell's equations 
for the modes Hz{x,y,t) — ip{x,y)e^^'^'^ reduce to a two- 
dimensional scalar wave equation [28:] 



(3) 



with frequency lo and the speed of light in vacuum c. The 
wave function ^ and its normal derivative times rt^^ are con- 
tinuous across the boundary of the cavity. At infinity, outgo- 
ing wave conditions are imposed. Even though the geome- 
try of the cavity is rather simple the wave equation cannot be 
solved analytically since the boundary conditions introduce 
diffraction at corners. We compute the modes numerically us- 
ing the boundary element method |29]. In order to apply this 
method each corner is replaced by a quarter of a circle with 
radius much smaller than the wavelength. We have carefully 
checked that the rounding does not influence the solutions in 
the studied frequency regime. 

Before we discuss the numerical solutions of the open cav- 
ity we will briefly consider the corresponding closed cavity 
with vanishing wave intensity along the boundary. We ex- 
pect that the solutions of the closed system approximate those 
modes in the open system which are confined by total internal 
reflection. The closed cavity is called an integrable billiard |l7J 
since the modes can be computed analytically 



^n^.ny{x,y) = sin(^^x) sin(^y) 



(4) 



if < a; < i? and < y < eR; otherwise ipn^,ny{x,y) = 
0. The positive integers Ur^, Uy count the number of maxima 
of \jp\'^ in X- and y-direction, respectively. The normalized 
frequency ft = ojR/c belonging to such a mode is given by 



As expected for an integrable billiard ff\, this system shows 
frequency crossings instead of avoided crossings when the as- 
pect ratio £ is varied. For example, for the modes (n^^, Uy) — 
(10,7) and (12,5) Eq. (|5]l yields the crossing point e = 
■y/6/11 « 0.739 and ft w 13.12 coiTesponding to a free- 
space wavelength of about 960nm. Figure |2] shows that this 
accidental degeneracy is lifted in the open cavity. The associ- 
ated ARC equals the case of the 2x2 matrix in Fig. [T] We 
therefore conclude that diffraction at corners in rectangular 
cavities leads to an external coupling of modes. 

It has been demonstrated in Ref. |30] that losses from a 
polygonal cavity due to diffraction at corners can be esti- 
mated by the boundary wave approach. Boundary waves 
travel along a flat interface between dielectric material and 
air. In the case of an infinitely extended interface, these waves 
are evanescent. In the case of a finite interface, however, these 
waves can leave the interface region at the corners. Follow- 
ing Ref. llsoll we have derived a formula describing the losses 
from a given mode (n^, Uy) in the rectangular cavity due to 
boundary waves 



lm{n) = 



2n 



2 

E 



sin 0A 



j=i 



(6) 



with a, 



n\ smf 



(5) 



J J -l/cos6j, tan 01 = eux/ny, 
and O2 — 7r/2 — 01. For a mode A with e = 0.72 and 
{nx,ny) — (10,7) we find Im(f2) w —0.0048 correspond- 
ing to a quahty factor of Q = Re(f2)/[2Im(17)] « 1380. 
For a mode B with e — 0.72 and {ux, riy) = (12, 5) we get 
Im(fi) « —0.0109 and Q « 600. The boundary wave ap- 
proach as developed in Ref. |30] can only compute the losses 
of individual modes, i.e. the diagonal elements of the ma- 
trix ([T]i. The off-diagonal part, i.e. the coupling of modes, 
cannot be determined within this approach. For a direct com- 
parison to the exact results in Fig. |2]it is therefore useful to 
consider the mean value of 17+ and ri_ since here the ARC 
contributions cancel; cf. Eq. Q. The result of this proce- 
dure is shown in the bottom panel of Fig. [3] It can be seen that 
the averaged boundary wave result overestimates the averaged 
lifetimes of the modes by just 20 percent. Hence, leakage due 
to boundary waves is the dominant decay channel. 

At the center of the ARC, e w 0.7453, in Fig.|2]a fast mode 
D with Im(ri) w -0.02 and a slow mode C with Im(r2) w 
-0.00028 is formed. The slow mode has Q « 23 200 which 
is a dramatic increase by more than one order of magnitude 
if compared to the "normal" quality factor. In this frequency 
regime the leakage due boundary waves limits the quality fac- 
tor to roughly 1900. This indicates that possibly all long-lived 
modes (modes with, say, Q > 4000) in this frequency regime 
are caused by ARCs. This conclusion is supported by exten- 
sive numerical studies on this system (not shown). 

The spatial patterns of modes A, B, E, and F in Fig. |4] ap- 
proximately match the solutions of the closed cavity in Eq. (01). 
Upon the avoided crossing the mode patterns exchange their 
character, i.e. mode B and E have roughly the same spatial 
profile but belong to different frequency branches; cf. Fig.|2] 
The same holds for mode A and F. The modes at the ARC, C 
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FIG. 2: An avoided resonance crossing in the rectangular microcav- 
ity. Plotted are the complex frequencies as function of the aspect 
ratio e. 
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FIG. 3: Mean value of real (top) and imaginary (bottom) part of 
frequencies (solid lines) near the avoided resonance crossing shown 
in Fig.|2] The dotted line is the scar approximation ([T}. The dashed 
line is the averaged result of the boundary wave approach l|6j. 



and D, correspond to symmetric and antisymmetric superposi- 
tions of the mode A and B (or E and F). Now, we can identify 
the physical mechanism behind the increased quahty factor of 
mode C: destructive interference reduces the light intensity at 
the corners and consequently the main decay channel, leak- 
age due to boundary waves, is strongly suppressed. A closer 
inspection of the mode C in Fig. |4] reveals that its intensity is 
concentrated along a diamond-shaped periodic ray. The long- 
lived mode formed in the ARC therefore resembles a scarred 
mode. In the case of mode D, constructive interference at the 
corners spoils the quality factor Mode D is localized along 
two symmetry-related rays connecting the corners of the cav- 
ity. Such kind of rays are called diffractive rays i^. 

The spatial pattern of mode C is not a special feature of 
the chosen boundary conditions but is also observed in square 
cavities of very different types: (i) square quantum dots with 
leads attached at the corner region UJ; (ii) square billiards 
with magnetic flux l IsTIl ; (iii) vertical-c avity surface-emitting 
lasers with square-shaped cross section Il3ll32ll . 

The resemblance of mode C with the diamond-shaped pe- 




••••••••••II 





)«ll«OII«l> 



lll l ll l l l ll 



•••••»•••••• 

••••loeiiiH 



■ 



FIG. 4: (Color online) Calculated near field intensity of modes with 
the same labels as in Fig.|2] Mode C shows localization along a peri- 
odic ray (dashed line), mode D is localized along two rays connecting 
the comers (dashed lines). 



riodic ray is apparent with the naked eye. In the following, 
it will be demonstrated that the relation is even deeper To 
do so, we estimate the frequency of the mode by using the 
localization along the ray. We stipulate that an integer num- 
ber m of wavelengths fits onto the periodic ray with length 
I = 2i?vl + e^. The calculation is straightforward giving 



= (m + /3) 



(7) 



J arctan aj being the total phase shift from 



with /3 = I 

the reflection at the dielectric boundary for TE polarization 
The quantities ai and a2 are the same as for Eq. @ but 
with tan 6*1 = and 62 = 7r/2 — 0i. The top panel of 
Fig.[3]demonstrates that the scar approximation with m = 15 
describes the mean behaviour of the modes involved in the 
ARC over a broad range of parameter values. The small fre- 
quency offset of about 0.086 can be traced back to the fact 
that the scar approximation assumes that the plane waves 
lying on the ray segments have no wave vector component 
in the transverse direction. However, the mode is restricted 
to an interval of length R (eR) in x-direction (y-direction). 
At least half a wavelength fits into these intervals for which 
the wave vector components give frequency contributions of 
approximately n/n and Tr/{en). Summing up the squares 
of the frequency contributions gives the correction Aft « 

(1 + l/e2) /2n2l7scai-- In the regime e = [0.72, 0.77] we 
get Aft « 0.096.. 0.102 which convincingly explains the dis- 
crepancy between fiscar and the exact value of fi. 

Coupling between two modes {nx,ny) and {mx,rny) oc- 
curs only for modes with the same symmetry with respect to 
the lines x = R/2 and y = eR/2. That implies that if 
is even (odd) rrix must be even (odd) too. The same holds 
for Uy and niy. Interestingly, this restriction ensures that for 
given allowed pair (n^^, Uy) and {rrix, my) cancellation at all 
corners is possible. 

We can create a variety of scarlike modes near ARCs. Con- 
sider a periodic ray bouncing q times at the horizontal lines 
and p times at the vertical lines. A straightforward analysis 



shows that q 
l)/2] where | 



2[(|n, 



-ma 



+ l)/2] andp = 2[(|nj 



denotes the integer part. Figure |3 a) depicts 
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FIG. 5: (Color online) a) Near field intensity of a long-lived mode re- 
sulting from an ARC of modes (10, 7) and (14, 5). The white dashed 
line is the periodic ray. b) and c) Long-lived TM mode formed at 
an ARC (marked by the dot) in an elliptical resonator with eccen- 
tricity e. d) and e) Long-lived TM mode at an ARC in a stadium 
resonator with horizontal lines of length L. 



an example with q — A and p = 2. This long-lived mode with 
« 16.306 - i0.00047 (Q « 17500) results from an ARC of 

modes (10, 7) and (14, 5) at e = 0.4954. 

The formation of long-lived, scarlike modes near ARCs 

with external coupling is not restricted to TE polarization nor 



to the rectangular geometry. Figures |3b)-(e) show examples 
with transverse magnetic (TM) polarization in an elliptical 
and a (semi-) stadium-shaped resonator with refractive index 
n = 3.3. The shape parameter is the eccentricity e and the 
length of the horizontal straight lines L, respectively. The 
scenario is as in the cases shown in Figs. [T] and |2l At the 
ARC a short- and a long-lived mode is formed; see Fig. ^c) 
and (e). The long-lived mode exhibits a localization along a 
periodic ray; see Fig. |3b) and (d). This localization grad- 
ually disappears when the shape parameter is detuned from 
resonance (not shown). Note that elliptical billiards do not 
show such scarring |33]. Let us mention that local maxima of 
quality factors as function of a shape parameter had already 
been exploited for minimizing losses from stadium-shaped 
cavities |12]. However, the case in Ref. [12] is not related 
to ARCs, but is an interference effect of unstable periodic 
rays Q. 

In summary, we demonstrated the formation of long-lived 
modes near avoided resonances crossings in optical microcav- 
ities. For a number of different types of cavities (rectangular, 
elliptical, and stadium-shaped) we observed strong localiza- 
tion of these modes, resembling scarred states. We expect that 
this finding is highly relevant for understanding the localiza- 
tion properties of long-lived states not only in optical systems 
but in various fields of research. 

We acknowledge helpful discussions with M. Hentschel, 
F. Anders T. Y. Kwon, and T. Gorin. 
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